HPLC and high-throughput sequencing revealed higher tea-leaves quality, soil fertility and microbial community diversity in ancient tea plantations: compared with modern tea plantations

Background Ancient tea plantations with an age over 100 years still reserved at Mengku Town in Lincang Region of Yunan Province, China. However, the characteristic of soil chemicophysical properties and microbial ecosystem in the ancient tea plantations and their correlation with tea-leaves chemical components remained unclear. Tea-leaves chemical components including free amino acids, phenolic compounds and purine alkaloids collected from modern and ancient tea plantations in five geographic sites (i.e. Bingdao, Baqishan, Banuo, Dongguo and Jiulong) were determined by high performance liquid chromatography (HPLC), while their soil microbial community structure was analyzed by high-throughput sequencing, respectively. Additionally, soil microbial quantity and chemicophysical properties including pH, cation exchange capacity (CEC), soil organic matter (SOM), soil organic carbon (SOC), total nitrogen (TN), total phosphorus (TP), total potassium (TK), alkali-hydrolyzable nitrogen (AN), available phosphorous (AP) and available potassium (AK) were determined in modern and ancient tea plantations. Results Tea-leaves chemical components, soil chemicophysical properties and microbial community structures including bacterial and fungal community abundance and diversity evaluated by Chao 1 and Shannon varied with geographic location and tea plantation type. Ancient tea plantations were observed to possess significantly (P < 0.05) higher free amino acids, gallic acid, caffeine and epigallocatechin (EGC) in tea-leaves, as well as soil fertility. The bacterial community structure kept stable, while fungal community abundance and diversity significantly (P < 0.05) increased in ancient tea plantation because of higher soil fertility and lower pH. The long-term plantation in natural cultivation way might significantly (P < 0.05) improve the abundances of Nitrospirota, Methylomirabilota, Ascomycota and Mortierellomycota phyla. Conclusions Due to the natural cultivation way, the ancient tea plantations still maintained relatively higher soil fertility and soil microbial ecosystem, which contributed to the sustainable development of tea-leaves with higher quality. Supplementary Information The online version contains supplementary material available at 10.1186/s12870-022-03633-6.


Background
Tea (Camellia sinensis L.) has become one of the most popular and widely consumed beverages in the world because of its multiple health benefits [1][2][3]. Southwest China including Yunnan, Guizhou and Sichuan Province has been regarded as the origin center, which is now widely cultivated all over the world for approximately 5000 years [4][5][6]. Until now, Lincang, Puer and Xishuangbanna Regions in Yunnan Province, China, still remained abundant ancient tea trees of wild, transitional and artificial culture types, respectively [7,8]. In particular, the fresh tea-leaves collected from the artificial culture trees with an age over 100 years in ancient tea plantations could be processed for Pu-erh tea, black tea and white tea [8].
Processing technology, tea cultivars, production places and harvest seasons significantly impacted phenolic compounds including catechins (a subgroup of flavan-3-ols), phenolic acids, flavonoids and anthocyanins, purine alkaloids, free amino acids like L-theathine, L-glutamate and L-arginine, and flavoalkaloids [9][10][11]. Additionally, compared with the ecological tea-leaves collected from modern tea plantations, white teas from the ancient tea plantations have been confirmed with higher ameliorative effect on kidney damage in diabetic mice via Sirtuins 1/AMP-activated protein kinase (SIRT1/AMPK) pathway, while Pu-erh teas from ancient tea plantations are discovered to possess higher contents of free amino acids, fatty acids, phenolic acids, nucleosides and nucleobases, but lower contents of flavonoids and caffeine congeners [12,13]. Generally, the fresh tea-leaves from ancient tea plantations showed relatively higher contents of tea polyphenols, total catechins, free amino acids, gallic acid, soluble sugars and caffeine, but significantly (P < 0.05) lower polyphenols/amino acids ratio and ester catechins content, especially epigallocatechin gallate (EGCG). Due to the lower yields, much more nutrients absorbed by deeper roots and greater chemical profile variations, the tea-leaves collected and processed from ancient tea plantations are widely recognized to have better quality and preserving value with a higher price [14].
The 15 rare earth elements contributed to the identification of ancient tea-leaves, which should be attributed to the soil nutrients and ecological system in ancient tea plantations [15]. As the supporter for survival and growth of tea trees, soil nutrients impacted the output and quality of tea-leaves, as well as microbial community composition [16]. In addition, the microbial including bacterial and fungal communities in rhizospheric soil play vital roles in ecological processes, such as soil organic matter (SOM) decomposition, soil nutrient cycling, and soil structure formation and stability [17,18]. Previous studies have indicated the variation of fungal community in rhizosphere soil with plantation age, and found the significant decrease of soil bacterial community in modern tea plantations after 23 years cultivation, which potentially caused the decrease of glomalin-related soil protein and soil organic carbon (SOC) after longterm cultivation [19,20]. Unfortunately, the distribution and characteristic of tea-leaves quality components, soil chemicophysical properties and microbial ecosystem in the ancient tea plantations remain unclear.
Thus, in this study, we evaluated the quality components of fresh tea-leaves, and soil chemicophysical properties and microbial community through high performance liquid chromatography (HPLC) and highthroughput sequencing in ancient and modern tea plantations from five geographical sites, which are all located at Mengku Town of Shuangjiang autonomous County in Lincang Region of Yunan Province, China. The following hypotheses were proposed: (1) tea-leaves quality components, soil chemical properties and microbial community would vary with the five geographical sites; (2) ancient tea plantation would differ from the modern tea plantation in soil chemical properties and microbial community, which potentially impacts quality composition of tea-leaves; and (3) soil environmental factors would impact bacterial and fungal community abundances and compositions due to the differences in tea plantation age and cultivation ways. Through comparisons, this study will help to elaborate the sustainable development of ancient tea plantations from soil fertility and microbial community structure.

Experimental sites
The study is performed at representative tea planting area that Mengku Town of Shuangjiang Autonomous County in Lincang Region of Yunan Province, which belongs to the subtropical mountain monsoon climate. The annual average temperature is 18 °C, and the annual average precipitation is from 1700 mm to 1900 mm. The soil type is Conclusions: Due to the natural cultivation way, the ancient tea plantations still maintained relatively higher soil fertility and soil microbial ecosystem, which contributed to the sustainable development of tea-leaves with higher quality. Keywords: Ancient tea plantation, Amino acids, Phenolic components, Soil fertility, Microbial community structure mainly laterite and yellow loam with the sand grains from granite sandy conglomerate and purple shale. Five sampling sites located at Mengku Town including Bingdao, Banuo, Baqishan, Dongguo and Jiulong all grow ancient tea plantations with an age over 100 years and modern tea plantations with an age no more than 60 years, respectively, which are present in Fig. 1. All tea plantations are grown a single cultivar of tea plant (Camellia sinesis var. assamica cv. Mengku Dayecha) identified by Yongmei Li and relevant experts, and cultivated in accordance with the requirement of ecological tea plantation that there were no farming, no fertilizers and no pesticides throughout the whole growing period.

Soil collection and treatment
The fallen leaves, weeds and stones on the surface of the soils from five subplots (0.5 × 0.5 m) in every plot (2 × 2 m) were clean up for the soil sampling collection. The soil samples were collected at the same cleared position. Sterile containers were used to store five soil samples at the 0-20 cm depth from five subplots in every plot, and these soil samples were blended into a composite soil sample. The field-moist clods (< 5 mm) underwent cool-drying at 4 °C until the gravimetric water concentration of about 80 g kg − 1 were separated through different sieve mesh including 1 mm, 0.25 mm and 0.15 mm for the determinations of soil physicochemical properties. The field-moist clods collected at 0-4 °C and − 86 °C in aseptic condition were prepared for the determination of microbiological quantification by culture medium and microbial community diversity by high-throughput sequencing, respectively.

Tea-leaves sample collection
The spring fresh tea-leaves with one bud and two leaves from different tea plantations were collected at the same period with reference to the method of Liu et al. (2020) [21]. The fresh tea-leaves were fixated by a 800 W microwave oven for 60-80 s, and then dried by sunlight until the moisture content below 6.0% by dry weight. The dried tea-leaves were maintained at − 20 °C for the chemical analysis, including tea polyphenols, free amino acids, 6 catechins, caffeine and gallic acid.

Determinations of main chemical components in tea-leaves
Moisture content and water extracts content were developed according to the national standard methods established by China National Institute of Standardization. Content of tea polyphenols was determined by a UV-1800PC ultraviolet-visible (UV-Vis) spectrophotometer (AOE Instruments, Shanghai, China) using Foline-Ciocalteu assay on 765 nm with gallic acid as the standard [22]. Total content of free amino acids was determined on 570 nm by the UV-Vis spectrophotometer using ninhydrin assay [22].
Six Catechins including catechin (C), epicatechin (EC), epigallocatechin (EGC), epigallocatechin gallate (ECG), gallocatechin gallate (GCG) and epigallocatechin gallate (EGCG), caffeine and gallic acid contents were determined by high-performance liquid chromatography B A (HPLC) through an Agilent 1200 series HPLC system (Agilent Technologies, Santa Clara, CA, USA) composed of diode array detector, infinity binary pump, integrated vacuum degasser, autosampler and thermostated column compartment and TSK gel ODS-80TM chromatogram column (250 mm × 4.6 mm, 5 μm; Tosoh Corporation, Tokyo, Japan) [23]. Solvent A (0.05 mol L − 1 phosphoric acid water solution and 5% acetonitrile) and solvent B (0.05 mol L − 1 phosphoric acid water solution and 80% acetonitrile) were prepared for HPLC separation. The gradient was programmed as follows:from 0 to 30 min, solvent A was decreased from 91 to 47%, solvent B was increased from 9 to 53% with a flow rate of 0.7 mL min − 1 . The column temperature and detection wavelength were set at 30 °C and 280 nm, respectively. Additionally, the injection volume was 5 μL. Quantitative analysis was carried out based on the calibration curves of 8 standards. Quantitative analysis was carried out based on the calibration curves of 8 standards. Each tea-leaves sample was determined with 5 replications.

Determinations of soil physicochemical properties
The soil pH value was measured by glass recombination electrode method in soil: water solution =1:2.5 (w/w). The cation exchange capacity (CEC) content was determined by ammonium acetate exchange extraction-distillation method. The soil organic carbon (SOC) content was determined by an UV-Vis spectrophotometer method through potassium dichromate with glucose as the standard. The soil organic matter (SOM) content was measured by potassium dichchromate capacity method-external heating method. The total nitrogen (TN) content was measured by Kjeldahl determination after digestion of concentrated sulfuric acid and hydrogen peroxide. The total phosphorous (TP) content was determined by molybdenum antimony anti colorimetric method after digestion of concentrated sulfuric acid and hydrogen peroxide. The total potassium (TK) content was determined by flame photometry after the digestion. The alkali-hydrolyzable nitrogen (AN) content was measured using alkalescent diffusion method of 0.5 mol L − 1 NaOH solution. The available phosphorous (AP) content was determined through phosphomolybdate blue spectrophotometry method after the extract of 0.5 mol L − 1 NaHCO 3 solution. The available potassium (AK) content was measured by flame photometry after the extract of ammonium acetate. Each soil sample was determined with 5 replications.

Determinations of microbial quantity in soil
The microbial quantity including fungi, bacteria and actinomycetes were determined by dilution coating method through the relevant culture mediums, and the colony forming units (CFU) were calculated by per gram of dried soil [24]. Beef extract peptone medium was carried out to measure bacterial quantity after 3 days of cultivation at 28 °C, while fungi number was measured through Martin's rose Bengal medium after 5 days of cultivation at 28 °C. And the actinomycetes in soil were determined by Gao's No.1 medium after 7 days of cultivation at 28 °C.
The fastp software was adopted for the quality control of original sequencing sequence. The FLASH software was carried out for DNA fragment assembly, through which sequences of < 50 bp and the reads with ambiguous bases, could be removed. The same operational taxonomic units (OTUs) were clustered by the sequences at the 97% identity threshold using UPARSE software. An 70% threshold was adopted to determine the taxonomic assignments of the OTUs through the Silva and Unite database. At OTUs level, the soil microbial community diversity, including the Chao 1 estimator and Shannon diversity were calculated using MOTHUR software [26].

Statistical analyses
The raw data of tea-leaves quality components, soil environmental factors and microbial community composition in modern and ancient tea plantations of five sampling sites are present in Additional file 1. The independent-samples T-test was carried out to compare tea-leaves chemical components, soil chemicophysical properties, and microbial community structure between ancient and modern tea plantations using SPSS 20.0 for Windows (IBM, Armonk, NY, USA). Principal coordinates analysis (PCoA) and non-metric multidimensional scaling (NMDS) were employed to evaluate the differences of soil microbial community composition according to Bray-Curtis distances between two types of tea plantations using Canoco 4.

Tea-leaves quality components variations with geographic site and tea plantation type
Except that the significant (P < 0.05) differences of chemical components were found in tea-leaves among five various geographic sites (Additional file 2: Table S1), the independent-samples T-test analysis indicated the definite differences of tea-leaves chemical components between modern and ancient tea plantations (Table 1), such as free amino acids and catechins, which revealed chemical components variations caused by geographical distribution and plantation type. In all sampling sites (Additional file 2: Table S1), tea-leaves collected from Bingdao possessed the highest levels of chemical components, including free amino acids, caffeine, gallic acid, GCG and ECG. Conversely, the tea-leaves from Baqishan had highest polyphenols/amino acids ratio value, and relatively higher contents of tea polyphenols, C and EGC, which was consistent with the report by   [28], who discovered the differences of metabolic components in raw Pu-erh teas from various mountains.
Among 16 chemical indicators, 11 chemical indicators showed significant (P < 0.05), highly significant (P < 0.01) or extremely significant (P < 0.001) differences, including tea polyphenols, free amino acids, caffeine, gallic acid, EGC, GCG and EGCG, particularly polyphenols/amino acids ratio and non-ester/ ester Table 1 Chemical indicators of tea-leaves in modern and ancient tea plantations, respectively Note: Polyphenols/Amino acids ratio = Tea polyphenols content/ free amino acids content; Total catechins were the summation of C, EC, EGC, ECG, GCG and EGCG contents; Non-ester catechins were the summation of C, EC and EGC contents; Ester catechins were the summation of ECG, GCG and EGCG contents; Non-ester/ Ester catechins ratio = Non-ester catechins content/Ester catechins content And all data were present by mean value ± standard deviation. The range was the mean value at 95% confidence coefficient. Significance difference levels: * P < 0.05; ** P < 0.01 and *** P < 0.001 determined by the independent-samples T-test using SPSS 20 catechins ratio with the extremely significant (P < 0.001) differences between modern and ancient tea plantations. Compared with modern tea plantations, ancient tea plantations were confirmed to possess higher contents of free amino acids, tea polyphenols, caffeine, gallic acid, EGC and total non-ester catechins, but had significantly (P < 0.05, P < 0.01 or P < 0.001) lower contents of GCG, EGCG and total ester catechins, which was consistent with the report by Cheng et al. (2018) [29] in EGCG content, while Ge et al. (2021) [13] discovered the catechins levels decreased after long-term plantation for about 1000 years. The extremely significantly (P < 0.001) lower polyphenols/amino acids ratio and higher non-ester/ ester catechins ratio in fresh tea-leaves from ancient tea plantations indicated the remarkable difference in chemical compositions distribution of free amino acids, catechins and phenolic compounds between modern and ancient tea plantations. Generally, the average contents of free amino acids, gallic acid and EGC in tea-leaves from ancient tea plantations were 1.22, 6.22 and 1.67 times higher than that from modern tea plantations, which should be attributed to the physiological metabolisms of tea plant, such as amino acid metabolism, caffeine metabolism and phenolic compounds metabolism associated with the soil environment and ecological system of tea plantation. Due to the nitrogen metabolism in tea plant, several new flavoalkaloids were found in green tea-leaves collected from the ancient tea plantation [29].

Soil chemicophysical properties and microbial quantity in various tea plantations
Soil chemicophysical properties and microbial quantity varied significantly (P < 0.05) with five geographical sites and two types of tea plantations (Table 2 and Additional file 2: Table S2). The highest levels of SOM, AN and AP were observed in Bingdao site, while Jiulong site possessed the relatively lower soil fertility (Additional file 2: Table S2). Conversely, the soil pH (ranging from 4.18 to 4.77) kept relative stable among five geographical sites. About microbial quantity, there were no significant (P ≥ 0.05) differences between two types of tea plantations on the whole. However, fungi population was extremely significantly (P < 0.001) increased in ancient tea plantations (11.25 ± 3.54 × 10 3 CFU g − 1 ) compared with the modern tea plantations (8.22 ± 1.63 × 10 3 CFU g − 1 ), which indicated that the long-term plantation age or natural cultivation way improved fungi community quantity in ancient tea plantations. However, bacteria and actinomycetes kept relatively stable with no significant (P ≥ 0.05) variations in modern and ancient tea plantations. Compared with the modern tea plantations, due to the micro-climate environment, low intensity production mode and less disturbance in cultivation management, the ancient tea plantations contained significantly (P < 0.05) higher soil fertility, including SOM, SOC, CEC, TN, TP, AN, AP and AK, but relative lower soil pH and TK. CEC (10.98-13.33 cmol kg − 1 , range across all the soil samples), SOC (47.89-54.73 g kg − 1 ), SOM Table 2 Soil physical properties and microbial population in modern and ancient tea plantations, respectively CEC Cation exchange capacity, SOC Soil organic carbon, SOM Soil organic matter, TN Total nitrogen, TP Total phosphorus, TK Total potassium, AN Alkali-hydrolyzable nitrogen, AP Available phosphorus, AK Available potassium All data were present by mean value ± SD. The range was the mean value at 95% confidence coefficient. Significance difference levels: * P < 0.05; ** P < 0.01 and *** P < 0.001 determined by the independent-samples T-test using SPSS 20 Table 2). The soil in modern tea plantations was high acidification after 60-years plantation [30]. Despite that the soil pH value (ranging from 4.38 to 4.55) in ancient tea plantations was significantly (P < 0.05) lower than that in modern tea plantations, it was still optimum pH condition for the growth of tea trees [31]. The higher CEC, SOM and SOC in ancient tea plantations delayed soil acidification trend through relevant microbial metabolic activity [32].
In spite of a long-term plantation age over 100 years, the ancient tea plantations still kept relatively higher soil microbial metabolic activities, particularly fungi communities with an extremely significant (P < 0.001) increase, which should be attributed to the natural cultivation way of ancient tea plantations [33]. Generally, the ancient tea plantations were cultivated close to tall trees and sheltered by virgin forests that contributed to the accumulation of SOM and SOC due to the higher litter derived from tea plant and other plants [34]. The higher nitrogen (N), phosphorus (P) and potassium (K) levels, particularly AN, AP and AK, ensure the sustainable development for the healthy growth of tea trees and the production of fresh tea-leaves in the ancient tea plantations. Except the TK with a relatively higher loss that was translated into AK under the effect of related microbial communities, the natural cultivation way was conducive to a growth over 100 years in ancient tea plantations because of the higher soil fertility. The higher AN, AP and AK improved theanine and flavonoid biosynthesis activities in tea plant physiology through glutamine synthetase, phenylpropanoid, flavonoid pathways [35,36], which might contribute to the accumulations of free amino acids, gallic acid, caffeine and EGC in tea-leaves. Additionally, the photosynthetic capacity potentially impacted catechins biosynthesis (such as EGC) in tea plant [37], which might cause the significant (P < 0.05) differences of EGC, EGCG and GCG between modern and ancient tea plantations.

Soil microbial community abundance and diversity in various tea plantations
A total of 6978 OUTs belonging to 30 phyla, 125 classes, 293 orders, 455 families, 856 genera and 1818 species were obtained by analyzing the bacterial communities, while a total of 5525 OUTs belonging to 20 phyla, 68 classes, 170 families, 853 genera and 1408 species were obtained by analyzing the fungal communities in all soil samples through the Illumina high-throughput sequencing, respectively. The values of bacterial coverage estimators ranged from 96.94 to 97.70%, and the values of fungal coverage estimators were all over 99.23% at 97% similarity cutoff, which suggested that the current numbers of the sequence reads were sufficient to reflect the soil bacterial and fungal community diversity (Additional file 2: Table S3).
Bacterial community abundance kept relatively stable among five geographic sites and two types of tea plantations ranging from 3375 OUTs to 4508 OUTs. Conversely, the independent-samples T-test analysis indicated that the fungal community abundance varied significantly (P < 0.05) with the geographic sites and tea plantation types. Generally, the ancient tea plantations involved a significant (P < 0.05) higher fungal community abundance. Therefore, plantation age enhanced the soil fungal community abundance in the natural cultivation way, while the bacterial community abundance kept stable with slight increases in the ancient tea plantations.
About α-diversity evaluated by Chao 1 and Shannon indices (Additional file 2: Table S3; Table S4), soil bacterial abundance kept stable with no significant (P ≥ 0.05) variations among five geographical sites and two types of tea plantations, except that the Chao 1 increased highly significantly (P < 0.01) in the ancient tea plantation (3215.36) of Jiulong site compared with the modern tea plantation (2754.08). However, Chao 1 and Shannon indices in soil fungal community diversity varied significantly (P < 0.05) with geographical sites and tea plantations type. Generally, α-diversity (i.e Chao 1 and Shannon indices) of fungal communities in ancient tea plantations were highly significantly (P < 0.01) or extremely significantly (P < 0.001) higher than that in modern tea plantations, respectively. Particularly, among five geographical sites, the ancient tea plantations located at Banuo, Dongguo and Jiulong sites showed significantly (P < 0.05) or highly significantly (P < 0.01) higher fungal community diversity than modern tea plantations. Compared with bacterial communities remaining stable along with the plantation age and cultivation ways, the fungal communities abundance and α-diversity significantly (P < 0.05) increased in the ancient tea plantations, which should be attributed to the long-term plantation age and natural cultivation way.
Previous studies confirmed that bacterial community diversity significantly (P < 0.05) decreased after longterm plantation over 23-years in modern tea plantations [19,30]. However, despite the long-term plantation over 100 years, the bacterial community abundance and diversity kept stable in ancient tea plantations because of the natural cultivation way that maintained relatively higher soil fertility. The significantly (P < 0.05) higher fungal abundance and diversity in ancient tea plantations was consistent with the report by Wu et al. (2020) [38], which confirmed that the long-term plantation might improve the fungal community in rhizosphere soil.

Bacterial community composition in modern and ancient tea plantations
Both PCoA ( Fig. 2A) and NMDS (Additional file 2: Fig.  S1) indicated the remarkable differences of bacterial community composition in soil aggregates between and ancient tea plantations. In the PCoA, the first two principal coordinates together represented 57.26, 70.54, 72.60, 73.03 and 77.31% of the variations by PC1 and PC2 in soil bacterial community composition from two types of tea plantations located at Bingdao, Baqishan, Banuo, Dongguo and Jiulong, respectively, which demonstrated that the bacterial community from same type of tea plantation in five various geographical sites was almost grounded together, while the bacterial community from the ancient tea plantation in the same geographical site were well separated from each other. Additionally, NMDS analysis revealed that the bacterial community composition between modern and ancient tea plantations in the same geographical site was significantly (P < 0.05) distinguished.
At the phylum level in soil bacterial community composition (Fig. 2B), Chloroflexi (24.94%), Proteobacteria (22.39%), Acidobacteriota (17.00%) and Actinobacteriota (13.15%) were regarded as the dominant groups in ancient tea plantations, which kept relatively stable in modern tea plantations with no significant (P ≥ 0.05) variations. Compared to the modern tea plantations, Proteobacteria and Actinobacteriota showed slight downward trend, while Chloroflexi and Acidobacteriota showed slight rising tendency in ancient tea plantations. Particularly, due to the long-term plantation and natural cultivation way, the relative abundances of Nitrospirota and Methylomirabilota significantly (P < 0.05) or highly significantly (P < 0.01) increased, but the relative abundances of Bacteroidota, Patescibacteria and Bdellovibrionota significantly (P < 0.05) or extremely significantly (P < 0.001) decreased in the ancient tea plantations (Additional file 2: Table S5).
The dominant phyla in bacterial community composition of tea plantations were consistent with the studies of  [30] discovered that Proteobacteria, Acidobacteriota, Actinobacteriota and Firmicutes were dominate phyla in soil of tea plantations. Although Actinobacteriota increased toward lower soil pH often [40], it kept stable in relative abundance in ancient tea plantations. As the main soil probiotics, the slight improvement of Bacillus genus enhanced the soil fertility and conduced to the sustainable development of ancient tea plantations.

Fungal community composition in modern and ancient tea plantations
Similar to bacterial community composition, fungal community composition varied significantly (P < 0.05) with geographical site and tea plantation types through PCoA (Fig. 3A) and NMDS analysis (Additional file 2: Fig. S2). In the PCoA, the first two principal coordinates together represented 55.90, 48.10, 65.79, 57.62 and 45.75% of the variations by PC1 and PC2 in soil fungal community composition between modern and ancient tea plantations located at Bingdao, Baqishan, Banuo, Dongguo and Jiulong, respectively, which indicated the remarkable differences in soil fungal community composition between modern and ancient tea plantations. The results of PCoA were consistent with NMDS analysis (Additional file 2: Fig. S2), demonstrated the extremely significantly (P < 0.001) distinction in fungal community composition between two types of tea plantations in the same geographical site.
At the phylum level in soil fungal community composition (Fig. 3B), Ascomycota (65.75%), Basidiomycota (19.02%), Mortierellomycota (5.08%), Rozellomycota (1.61%), Glomeromycota (0.39%) and Unclassified-k-Fungi (7.57%) were the dominant groups in ancient tea plantations. Except for Jiulong, the relative abundance of Ascomycota in the ancient tea plantations of other four sites was significantly (P < 0.05) higher than that in modern tea plantations, which occupied from 51.50 to 70.28%. The soil fungal community composition were significant (P < 0.05) differences across five sampling sites. For example, the relative abundance of Rozellomycota in modern tea plantation of Dongguo site reached 20.77%, but its relative abundance in other soil aggregates was relatively low. Compared to modern tea plantations, Unclassified-k-Fungi and Rozellomycota significantly (P < 0.05) decreased, while Ascomycota and Mortierellomycota significantly (P < 0.05) increased in the ancient tea plantations (Additional file 2: Table S7).
The dominant phyla in fungal community composition were consistent with the report by Wu et al. (2020) [38]. The abundances of Ascomycota and Mortierellomycota phyla, as well as 9 genera (i.e. Saitozyma, Metarhizium, Mortierella, f__Clavariaceae, Pseudogymnoascus, o__Helotiales, Leohumicola, Tolypocladium and Trichoglossum) significantly (P < 0.05) increased with the plantation age in natural cultivation way, which should be attributed to the higher fertility in the  (PCoA, A) showed the remarkable difference of soil fungal community composition between modern and ancient tea plantations in five various sampling sites, and their relative abundance differences in the top 12 phyla (B) through independent-samples T-test, respectively ancient tea plantations. The entomopathogenic fungi stay at relatively lower level in ancient tea plantations. For instance, the abundance of Fusarium genus showed a slight decrease from 2.86 to 1.82% in the ancient tea plantations.

Potential impact of soil physical properties in tea plantations
The bivariate correlation analysis was carried out to explore the potential impact of soil physical properties on tea-leaves quality components and microbial community in tea plantations through the relevant Spearman correlation coefficient (Additional file 2: Table S9; Table S10). The soil TN, TP, AN, and AP levels might improve the contents of free amino acids, gallic acid and EGC in tea-leaves with the significantly (P < 0.05) or highly significantly (P < 0.01) positive correlations (r-value ranging from 0.355 to 0.624), but decreased EGCG content with the highly significantly (P < 0.01) negative correlations (r-value = − 0.451, − 0.471 and − 0.611, respectively). Additionally, TK and AK impacted potentially the levels of phenolic components including gallic acid, C, EC, EGC, ECG, GCG and EGCG, due to the significantly (P < 0.05) or highly significantly (P < 0.01) negative correlations (Additional file 2: Table S9). These results indicated that soil environmental factors, such as CEC, SOC, SOM, nitrogen (N), phosphorus (P) and potassium (K), showed significant (P < 0.05) associations with tealeaves quality components, particularly free amino acids and phenolic components. Generally, SOC, SOM, N and P might contribute to the accumulations of free amino acids, gallic acid and EGC, while the K might impact the biosynthesis of phenolic compounds, including gallic acid, C, EC, EGC, ECG, GCG and EGCG, which generated the chemical differences in tea-leaves between modern and ancient tea plantations.
The soil pH was main environmental factor impacting the microbial abundance and diversity. In the modern tea plantations, the bacterial abundance and diversity decreased significantly (P < 0.05) after the long-term plantation due to the continuous acidification in the soil [19,28]. However, despite the reduction of bacterial abundance with the lower soil pH, the relatively lower soil pH improved the abundances of actinomycetes and fungi in ancient tea plantations. Additionally, the microbial diversity including bacterial and fungal diversity was not sensitive to soil pH (Additional file 2: Table S10). Differed form soil pH, the soil fertility, such as SOC, SOM and TP, might improve fungi abundance determined by dilution coating method. Additionally, the soil P level showed the significantly (P < 0.05) or highly significantly (P < 0.01) positive correlations with the abundance and diversity of fungi community, and actinomycetes richness, while AK level might improve bacterial community abundance. The N level in soil had limited influence on soil microbial community abundance and diversity with no significant (P ≥ 0.05) correlations to bacterial and fungal communities.

Interactions between soil physical properties and microbial community composition
The Spearman correlation coefficient determined by the bivariate correlation analysis corresponded highly to the RDA. In this study, the top 22 bacterial phyla (Fig. 4A) and 12 fungal phyla in microbial communities were carried out as the response variables through the RDA to elaborate the impact of soil environmental factors on microbial community composition in tea plantations, respectively. The RDA model accounted for 93.56 and 86.87% of the total variation of soil bacterial and fungal community composition in tea plantations, respectively. AK, AN, TP, CEC and TK were main environmental factors affecting bacterial community composition in turn (Fig. 4A). CEC, TK and soil pH showed significantly (P < 0.05) positive correlations with the relative abundances of Proteobacteria and Acidobacteriota, while the positive correlations of the relative abundances of Actinobacteria and Chloroflexi to AK, AN, and AP were found in the bacterial community composition (Fig. 4A). Therefore, CEC contributed to the improvement of f__Xanthobacteraceae, while TK and soil pH impacted Sphingomonas and Burkholderia-Caballeronia-Paraburkholderia genera in ancient tea plantations. Additionally, the soil AN, AP and AK showed positive correlations with the abundances of c__TK10, HSB_OF53-F07, Mycobacterium and Kitasatospora in ancient tea plantations.
However, for the fungal community composition in tea plantations, TP, AN, AP, pH and AK were principal influence factors. The abundance of Ascomycota phylum in fungal community composition showed significantly (P < 0.05) positive correlations to soil chemical properties, such as AN, AP, AK, SOC, SOM and TP, but with negative correlation to soil pH, which indicated that the soil chemical properties improved the abundance of Ascomycota phylum including Metarhizium, Pseudogymnoascus, o__Helotiales, Leohumicola, Tolypocladium and Trichoglossum genera toward to the lower pH in the ancient tea plantations. Soil pH and TN showed definitely positive corrections to the abundances of Unclassified_k__Fungi, Basidiomycota and Rozellomycota, which revealed that these fungal phyla were sensitive to the soil pH with significant (P < 0.05) decreases in ancient tea plantations (Fig. 4B). Generally, the significant (P < 0.05) increases of Ascomycota phylum as well as Metarhizium, Pseudogymnoascus, o__Helotiales, Leohumicola, Tolypocladium and Trichoglossum genera in the abundances were caused by the higher soil fertility, and lower pH value contributed to the stabilization of soil microbial diversity and abundance in ancient tea plantations.
Conclusions Due to the excellent ecological environment, natural cultivation way and complete community structure, the ancient tea plantations were observed to possess significantly (P < 0.05) higher free amino acids, caffeine, gallic acid and EGC contents, but lower GCG and EGCG contents in tea-leaves, as well as higher soil fertility than modern tea plantations, including CEC, SOC, SOM, TN, TP, AN, AP and AK. The long-term plantation in natural cultivation way significantly (P < 0.05) improved the abundances of Nitrospirota and Methylomirabilota phyla as well as 5 genera in bacterial community composition, and the abundances of Ascomycota and Mortierellomycota phyla as well as 9 genera in fungal community composition. The bivariate correlation analysis and RDA revealed the impact of soil physical properties on tea-leaves chemical components and microbial community composition in tea plantations. The significantly (P < 0.05) higher soil fertility promised the stabilization of bacterial abundance and diversity, and improved the abundance and diversity of fungal community. Generally, This study revealed the characteristics of tea-leaves chemical components, soil chemicophysical properties and microbial community structure in the ancient tea plantations, which advances the knowledge about the sustainable development of ancient tea plantations in soil ecosystem.